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Effect of geometric and electronic structures 
on the finite temperature behavior of Nasg, Nasy, and Nass clusters 

Mal-Soon Le^EI and D. G. KanhertJl 

Centre for Modeling and Simulation, and Department of Physics, 

University of Pune, Ganeshkhind, Pune ^11 007, India. 
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An analysis of the evolutionary trends in the ground state geometries of Nass to Na62 reveals Nasg, 
an electronic closed-shell system, shows namely an electronically driven spherical shape leading to 
a disordered but compact structure. This structural change induces a strong connectivity of short 
bonds among the surface atoms as well as between core and surface atoms with inhomogeneous 
strength in the ground state geometry, which affects its finite-temperature behavior. By employing 
ab initio density-functional molecular dynamics, we show that this leads to two distinct features 
in specific heat curve compared to that of Nass: (1) The peak is shifted by about 100 K higher 
in temperature. (2) The transition region becomes much broader than Nass. The inhomogeneous 
distribution of bond strengths results in a broad melting transition and the strongly connected 
network of short bonds leads to the highest melting temperature of 375 K reported among the 
sodium clusters. Nasr, which has one electron less than Nass, also possesses stronger short- bond 
network compared with Nass, resulting in higher melting temperature (350 K) than observed in 
Nass. Thus, we conclude that when a cluster has nearly closed shell structure not only geometrically 
but also electronically, it show a high melting temperature. Our calculations clearly bring out the 
size-sensitive nature of the specific heat curve in sodium clusters. 

PACS numbers: 31.15.Ew,36.40.Cg,36.40.Ei,64.70.Dv,71.15.Mb 



I. INTRODUCTION 

The finite-temperature behavior of sodium clusters has 
been attracted much attention since the pioneering ex- 
perimental work by Haberland and co-workerssi^ These 
experiments reported melting temperatures of sodium 
clusters in the size range of 55 to 350 atoms. A great deal 
of effort has been spent to understand and explain the 
main puzzling feature, namely the irregular variation in 
the observed melting points. Equally puzzling is the ob- 
servation that Na57 and Nai42 , which is neither geometric 
shell-closing systems nor electronic shell-closing systems, 
show higher melting temperatures than vicinity systems 
such as Nass and Nai47 of geometric shell-closing ones 
or Nai38 of electronic shell-closing one, which is not ex- 
plained so far. Although it is expected that there is an 
intricate interplay between the electronic and geometric 
structures so far as finite temperature behavior is con- 
cerned, the consensus seems to be that the geometry 
dominates the melting characteristics and the effects of 
electronic structure are secondary:?'-?- This has been a sur- 
prise since the stability of these clusters is almost entirely 
dictated by the electronic structure. According to our 
present calculations, it reveals that when a clus- 
ter with nearly electronic closed shell structure 
has nearly icosahedral ground state geometry, it 
shows a high melting temperature. 

Recently a very different aspect of the finite- 
temperature behavior has been revealed in the experi- 
mental measurements of heat capacities of gallium and 
aluminum clusters;* namely the size-sensitive nature of 
their shapes. A definitive correlation between the shape 
of specific heat curve and the nature of the ground-state 



geometry in gallium cluster has also been established by 
our group. ^"^ Surprisingly, in spite of extensive work re- 
ported on the bench mark system of sodium clusters, 
there is no report of such a size sensitivity. 

In this paper, we show the effect of electronic struc- 
ture on the melting of sodium cluster by comparison to 
a geometric shell-closing system. It is interesting to note 
that there are very few cases in sodium clusters where an 
electronic shell-closing system (N=8, 20, 40, 58, 138, ...) 
has a similar size as a geometric shell-closing one (N=55, 
147, 309, ...). For instance, a pair of clusters N=55 and 
N=58 differs by three atoms only. Another example is a 
pair {N=138 and 147} having nine-atom difference. The 
next such occurrence is at N=309 and N=338, which is 
not close enough. Therefore, the expected effects due 
to both the geometric and electronic magic numbers are 
likely to be seen prominently in {55, 58} range. The 
electronically closed-shell structure of Nasg and nearly 
electronic shell closing system of Nasy have slightly dis- 
torted icosahedral structures as their ground state ge- 
ometries compared to the geometric shell-closing cluster 
of Nass. Indeed, as we shall see, Nass and Nasr, have 
significantly different finite-temperature characteristics 
as compared to those of geometrically closed-shell clus- 
ter, Nass. Firstly, the shape of the specific heat curve 
is much broader than that of Nass- Secondly, it shows 
a peak (rather broad) at the temperature approximately 
375 K for Nass and 350 K for Nas7. We note that the 
observed melting temperature of Nas7 by Haberland and 
co-workers '•. is about 325 K which differs less than 8 % 
from our result. Specially, in case of Nass this is the first 
observation that a melting temperature of sodium clus- 
ter can be closed to that of the bulk. Interestingly, it is 



observed that Na58 shows high abundance in mass spec- 
tra, indicating its high stabihty,^ The paper is organized 
as follows: in Sec. Ej we briefly mention computational 
details used. In Sec. lIIII we discuss our results, and sum- 
marize the results in Sec. IIVI 



II. COMPUTATIONAL DETAILS 

We have carried out Born~Oppenheimer molecular- 
dynamics simulations using ultrasoft pseudopotentials 
within the local-density approximation (LDA)i& The re- 
liability of our calculation can be judged from the fact 
that our earlier calculations based on density-functional 
theory (DFT) have successfully reproduced the melting 
temperature of Na„(n=55, 92, and 142)ii In order to get 
an insight into the evolutionary pattern of the geome- 
tries, we have obtained equilibrium structures for Na„ 
(n=55 to 62). Our thorough search of the lowest-energy 
structure is done by obtaining at least 180 distinguish- 
able equilibrium conflgurations for each of the clusters 
by using a combination of a basin-hopping algorithm 
and density-functional methods. With obtained ground- 
state geometry of Nasg and Nasy, we have carried out 
extensive ab initio constant-temperature simulations us- 
ing a Nose-Hoover thermostat to compare with geomet- 
rically closed shell system, Nass. These simulations have 
been carried out at 16 temperatures for Nasg and at 12 
temperatures for Nasy in the range of 80K< T<500K 
for the period of at least 210 ps (240 ps near the melting 
temperature), total simulation times of 2.7 ns and 3.5 ns, 
respectively. Our cell size used is 24x24x24 A'^, with the 
energy cutoff of 3.6 Ry whose reliability has been exam- 
ined in our previous work.^ With obtained ground-state 
geometry of Nasg and Nasy, we have carried out extensive 
ab initio constant -temperature simulations using a Nose- 
Hoover thermostat to compare with geometrically closed 
shell system, Nass. These simulations have been carried 
out at 12 temperatures for Nasy and at 16 temperatures 
for Nags in the range of 80K< T<500K for the period 
of at least 210 ps (240 ps near the melting temperature), 
total simulation times of 2.7 ns and 3.5 ns, respectively. 
Our cell size used is 24x24x24 A'^, with the energy cutoff 
of 3.6 Ry whose reliability has been examined in our pre- 
vious workji We use the multiple-histogram technique to 
calculate specific heat. More details of the methods can 
be found in Ref. la- For the analysis, we have taken the 
last 165 ps data from each temperature, leaving at least 
first 45 ps for thermalization. 



III. RESULTS AND DISCUSSION 

The lowest-energy structures of Na„ (n=55-62) are 
shown in Fig. Q] There are some striking features evi- 
dent in their evolutionary pattern. The ground-state ge- 
ometry (GS) of Nass is the well-known icosahedron^. A 
single extra atom added to Nags is accommodated on the 
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FIG. 1: The ground-state geometry of Na„ (n=55-62). In 
the figure of Nasg , "A" indicates the region where three extra 
atoms are accommodated compared to Nass, and "B" indi- 
cates the region where the geometry of Nass is retained. 
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FIG. 2: The distance from the center of mass of Nass and 
Nasg in comparison with that of Nass . In Nasg the maximum 
distance is almost the same as seen in Nass, while that in 
Nas9 is changed significantly. 



surface by a minor adjustment of surface atoms. When 
the second atom is added, it is energetically more fa- 
vorable to retain the icosahedral core with two atoms 
capping it. Very interestingly, when three atoms are 
added (n=58), instead of the pattern of capping con- 
tinuing, all three atoms are accommodated on the sur- 
face of icosahedron, making the structure again nearly 
spherical without changing the size of Nass significantly. 
This can be verified by examining the distance from the 
center of mass (COM). In Fig. [3 we show the distance 
from COM in Nass and Nasg by comparison with those 
of Nass. It is clearly seen that in comparison with the 
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FIG. 3: The comparison of (a) the deformation parameter, 
Edef, (b) the eigenvalue spectra of the ground-state geome- 
tries as a function of cluster size. 



maximuni distance from COM in Nass, that in Nasg is 
nearly remained the same, while in Nasg it increases con- 
siderably. The pattern of growth from Nasg to Na62 
changes back to the capping mode as can be seen in 
Fig. n This peculiar shape transformation observed in 
Nass can be examined by plotting deformation param- 
eter Edef- £def is defined as Edef = '^Qil{Qi + Qa) , 
where Q\ > Q2 > Q3 are eigenvalues of the quadrupole 
tensor Qij — J^i RuRij with Rn being «*^ coordinate 
of ion / relative to the center of mass of the cluster. For 
a spherical shape (Qi = Q2 = Qz) £def is 1.0, while 
Edef > 1-0 indicates a deformation. It can be seen in Fig. 
Ola) that the addition of three atoms over Nass changes 
the shape to nearly spherical. Interestingly, this differ- 
ence is reflected in their eigenvalue spectrum as shown in 
Fig. |2Ib) . A peculiarity of Nass is that the structure is 
not highly symmetric, rather in the sense of amorphous. 
However, it follows jellium-like pattern very closely in 
contrast to Nasr, Nasg to Na62 where additional states 
appear in the gaps due to their disordered structures. 
Such an effect is absent in Nass- In addition an elec- 
tronic shell-closing system of Nass shows the highest gap 
of the highest occupied molecular orbital (HOMO) and 
the lowest unoccupied molecular orbital (LUMO) among 
present studied systems. Thus, an electronically closed- 
shell nature of Nass results into a spherical charge density 
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FIG. 4: The histogram of bond lengths for the ground state 
geometry of Nasg, Nasr, and Nass, with less than that of the 
bulk (3. 71 A). 



distribution which drives the geometry towards a spher- 
ical shape. Our calculations bring out the fact that this 
change of shape induces the GS geometry to be compact 
with a significant structural disorder in Nass. This is an 
example of an electronically driven shape change. 

We find that to understand the melting behavior of 
small cluster, which we shall see later, it is important to 
study not only how many short bonds are there in the GS 
geometry of the cluster, but also how these short bonds 
are connected each other (we call this ^^ connectiviti/^ ) . 
Thus, we have examined the bond lengths in Nass, Nasy, 
and Nass as well as their connectivity. Fig. 0] shows the 
number of bonds having bond lengths less than 3.71 A, 
the bulk bond length. Nass has 21 bonds shorter than 
the shortest bond in Nass , which are located in the region 
A, shown in Fig. ^ giving rise to an island of relatively 
strongly bonded atoms as compared to Nass- In Nas7 
there arc 16 bonds which are shorter than the shortest 
bond in Nass- Fig. [S] shows how the nature of connec- 
tivity in these clusters with different bond lengths. In 
Figs. EIa)-|5Ic), connectivity of short bonds are shown 
with bond lengths less than 3.45 A, and Figs. intd)-|SIf) 
for less than 3.55 A. It can be seen that Nass has fewer 
short bonds (12 bonds) than Nass (41 bonds) and Nas7 
(38 bonds). The majority of these strong bonded atoms 
in Nass form a connected island in the region A shown in 
Fig.^ while all short bonds in Nass are radial. It is also 
evident from Fig. \^d) that for bond lengths less than 
3.55 A Nass shows not only that core atoms are strongly 
bonded to surface but that even surface atoms are bonded 
each other. In contrast to this, there is a connectivity 
only between core (first shell) and surface atoms but not 
among surface atoms in Nass (Fig- 13^))- Nas7 shows 
stronger and inhomogeneous connectivity than those seen 
in Nass but weaker than those seen in Nass in both bond 
length regimes. Thus, we conclude that to accommodate 
three extra atoms in Nass first-shell and surface distance 
is reduced, resulting in a strong network extending over 
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FIG. 6: The isovalued surfaces of the electron locahzation 
function at xelf=0.79 for Nasg (a) in region A and (b) in 
region B. The similar pattern of (b) is also seen at xelf=0.76 
in Nass. 



FIG. 5: The short-bond connectivity with bond length 
less than 3.45A in (a) Nass (b) Nasy (c) Nass, showing 
inhomogeneous distribution of bond strength in Nass and 
Nas7. The short-bond connectivity with bond length less 
than 3.55A in (d) Nas8 (e) Nasr (f) Nass, showing Nasg has 
strongest connectivity through entire system among three 
clusters. 



entire cluster with inhomogeneous strength. To compare 
the GS structures of Nass and Nass, We also analyze the 
distribution of short bond lengths in Nass and Nass by 
comparing the distance from different shells. It turns out 
that while Nass has shorter bond length between center 
atom and first-shell, the average bond length between 
first-shell and outer-shell is shorter in Nass- In a re- 
cent work by Aguado and Lopez;^ they have attributed 
the higher melting temperature, seen in experiments of 
sodium clusterSji^ to the existence of shorter bonds be- 
tween the surface and first-shell atoms. Our calculations 
are consistent with this observation. 

We have examined the electron localization function 
(ELF). The ELFia is defined as 
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Here p = p{r) is the valence electron density and ipi^s are 
the KS orbitals. The Xelf is 1-0 for perfect localized 
function and 0.5 for plain waves. In Fig. El we show the 
ELF isosurface of Nass for isovalue of X-E-L-F=0.79 in two 



different regions. Figs.EJa) andj^b) depicts the ELF at 
0.79 in regions A and B, respectively. The contrast is ev- 
ident. In region A which consists of atoms connected by 
shorter bonds, the ELF isosurface forms a connected net- 
work showing the existence of a strongly bonded region 
containing at least 20 atoms. In contrast, in region B 
such a network is absent. In fact, the region B gets con- 
nected at a lower isovalue of 0.74. The same analysis has 
been carried out for Nsv and Nass. Nasr shows similar 
characteristics compared to those seen in Nass ■ A region 
near the capping atoms (left-hand side in Fig.|3Ib) starts 
connecting at xelf=^-^^, whereas a region away from 
them is connected at 0.74. The connection of isosurface 
is established over the entire cluster of Nass at the iso- 
value of 0.72. Thus, our ELF analysis clearly brings out 
the existence of strongly bonded island of atoms and the 
inhomogeneous distribution of bond strength in Nass and 
NasT. 

Our study reveals two unique features in the ground 
state geometries of Nass and Nas?. Firstly, the (nearly) 
electronic shell closing system induces shortening of 
bonds, resulting in a strong connectivity of short bonds 
among the surface atoms as well as between first-shell 
and surface atoms compared to only first-shell to core 
connectivity established in Nass. Secondly, its ground- 
state geometry is considerably disordered, resulting in in- 
homogeneous distribution of bond strength as compared 
to highly symmetric Nass. Hear, disordering in these 
cases is related to the absence of spherical symmetry. We 
argue that these two effects will be manifested in the spe- 
cific heat differently. The existence of the well-connected 
short-bonds are expected to raise the melting tempera- 
ture as compared to that of Nass • The effect of geometric 
disorder is to broaden the specific heat. 

We show the calculated specific-heat curves for Nass 
and Nas7 along with those of Nassr^ in Fig- El The most 
symmetric cluster Nass shows a sharp melting transi- 
tion at 280 K, while the highest melting temperature is 
seen in Nass. Thus, the addition of two or three atoms 
changes the specific heat curve drastically. Our calcu- 
lations clearly demonstrate a strong correlation between 
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FIG. 7; The normalized specific lieat as a function of temper- 
ature. Co ~ {SN — 9/2)kB is the zero-temperature classical 
limit of the rotational plus vibrational canonical specific heat. 



the nature of order or lack of order in the GS geome- 
try and their finite temperature behavior, and between 
the connectivity of short bonds and its melting temper- 
ature. This explanation can be extended to the finite 
temperature behavior in Naso (figure not shown). Naso 
has a disordered GS geometry but does not possess a 
strongly connected network.'''-^ The specific-heat curve 
of Naso, indeed, shows a very broad melting transition 
with low melting temperature (225 K). The results of 
the density-functional simulations carried by Rytkonen 
et alm^ as well as Lee et ali^ are also consistent with the 
features that emerge out of our present calculation. An 
electronically closed-shell system Na4o has a disturbed 
spherical GS geometry with a well-connected network of 
short bonds from core to surface atoms, which is stronger 
than that of Nass. Its specific heat is broader but melt- 
ing point is higher than that of Nass (figure not shown) . 
Thus, it also brings out the size-sensitive nature of spe- 
cific heat curve in sodium clusters, not reported so far ei- 
ther experimentally or computationally. It may be noted 
that extreme size sensitivity of this kind has been ob- 
served experimentally in Ga and Al clustersi^ We have 
also calculated the latent heats per tom by using caloric 
curves. They are estimated as 0.02 eV/atom for Nass, 
0.016 eV for Nasy, and 0.014 eV/atom for Nags- The 
latent heats observed by experiment^ are 0.008 eV/atom 




100 300 100 300 

Temperature (K) 



500 



FIG. 8: The comparison of the root-mean-square bond 
length fluctuation (Srms) of individual atoms for Nass and 
Nass. 



for Na57 and 0.015 eV/atom for Nass. The deviation 
observed in Nasy may be due to its very broad melting 
transition. 

Due to the disordered nature of the GS geometry in 
Nass, we expect different atoms to move with various 
amplitudes at a given temperature. To see this, we 
have calculated the Lindemann criteria, i.e., the root- 
mean-square bond length fluctuation {6rms) for individ- 
ual atoms of Nass, where Smisii) for atom i is defined 
as 
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where Rij is the distance between i and j ions with i ^j, 
N = n — 1 with n be number of atoms in the cluster, 
and {...)t denotes a time average over the entire trajec- 
tory. Indeed, 6rms(i) in Nass has a broad range of values 
at a given temperature showing a different response for 
each atom, while for Nass (highly symmetric structure) 
they are collective, leading to a sharp transition region. 
Further, it saturates at nearly 375 K in Nass compared 
to 280 K in Nass. Interestingly, the center atom in Nass 
(the bottom most line) dose not develop its melting be- 
havior till 300 K. The disordered nature of the geometry 
leading to inhomogeneous distribution of bond strength 
results in distinctly different response of different atoms 
upon heating. We show the mean-square-displacement 
(MSD) of individual atoms at 240 K in Fig. El Evidently 
the atoms in these two clusters respond differently when 
heated. While Nass shows a clear solid-like behavior, a 
significant number of atoms in Nass show diffusive mo- 
tion. We also observed the similar characteristics in Nasv, 
seen in Srms and MSD of Nass • 

As shown in Fig. ^| where we calculate the average 
eigenvalues carried over the entire simulation runs of 
Nass, Nas7, and Nass, we do not find that liquid-like 
and solid-like behavior of these clusters are electroni- 
cally similar. Instead, the HOMO-LUMO gaps in all the 
cases are closed. According to the calculations of Sdef 
as a function of temperature, after melting their geome- 
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FIG. 9: The comparison of the mean square displacement 
(MSD) of individual atoms at 240 K for Nasg and Nass. 
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tries become elongated in these clusters. This can explain 
the close of HOMO-LUMO gap. Thus, the present ex- 
tensive simulations clearly bring out the importance and 
significance of the electronic structure. Interestingly, the 
change of eigenvalue spectrum and Sdef as a function 
of temperature in Nasgand Na57 becomes similar after 
180 K. Even though the GS geometry of Nasv is not 
spherical, its structure changes towards spherical shape 
upon heating. Thus, when a cluster has both character- 
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FIG. 11: The average deformation parameter of Nass, Nas?, 
and Nasg as a function of temperature. 



istics of (nearly) closed shell not only geometrically but 
also electronically, it has high melting temperature com- 
pared to those having either geometric close shell or elec- 
tronic closed shell. This may explain the reason Nai42 
has high melting temperature than Naiag, an electronic 
closed-shell cluster and Nai47, a geometric closed-shell 
one. 



IV. SUMMARY 

To summarize, the present work brings out the effect 
of the electronic as well as the geometric structures on 
the melting of Nasg and Nasy clusters compared with 
that of Na55. The electronic shell-closing nature of Nasg 
drives the GS geometry to be spherical, which leads to 
a compact and disordered structure. As consequence, in 
Na5g the first -shell and surface as well as atoms on the 
surface are well connected with short bonds. This leads 
to the high melting temperature of nearly bulk melting 
temperature. The disordered nature of system is respon- 
sible for rather broad specific heat curve. Another ex- 
ample of strong connectivity is seen in the GS geome- 
try of Na57, nearly electronic shell closed cluster, which 
also shows high melting temperature than geometrically 
closed shell system of Nass. The strong correlation be- 
tween connectivity of bonds and melting temperature is 
also seen in the case of Gai2 and Al clusters.^'* Thus, we 
conclude that electronic structure affect its melting be- 
havior strongly. We believe that it is possible to verify 
the prediction of present work experimentally with the 
calorimetry method. 
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